Search for: All records

Abstract A consequence of a nonzero occupation fraction of massive black holes (MBHs) in dwarf galaxies is that these MBHs can become residents of larger galaxy halos via hierarchical merging and tidal stripping. Depending on the parameters of their orbits and original hosts, some of these MBHs will merge with the central supermassive black hole in the larger galaxy. We examine four cosmological zoom-in simulations of Milky Way-like galaxies to study the demographics of the black hole mergers that originate from dwarf galaxies. Approximately half of these mergers have mass ratios less than 0.04, which we categorize as intermediate mass ratio inspirals, or IMRIs. Inspiral durations range from 0.5–8 Gyr, depending on the compactness of the dwarf galaxy. Approximately half of the inspirals may become more circular with time, while the eccentricity of the remainder does not evolve. Overall, IMRIs in Milky Way-like galaxies are a significant class of black hole mergers that can be detected by LISA, and must be prioritized for waveform modeling. 

Abstract The interaction between supermassive black hole (SMBH) feedback and the circumgalactic medium (CGM) continues to be an open question in galaxy evolution. In our study, we use smoothed particle hydrodynamics simulations to explore the impact of SMBH feedback on galactic metal retention and the motion of metals and gas into and through the CGM of L_*galaxies. We examine 140 galaxies from the 25 Mpc cosmological volume Romulus25, with stellar masses between log(M_*/M_⊙) = 9.5–11.5. We measure the fraction of metals remaining in the interstellar medium (ISM) and CGM of each galaxy and calculate the expected mass of each SMBH based on theM_BH–σrelation (Kormendy & Ho 2013). The deviation of each SMBH from its expected mass, ΔM_BH, is compared to the potential of its host viaσ. We find that SMBHs with accreted mass aboveM_BH–σare more effective at removing metals from the ISM than undermassive SMBHs in star-forming galaxies. Overall, overmassive SMBHs suppress the total star formation of their host galaxies and more effectively move metals from the ISM into the CGM. However, we see little to no evacuation of gas from the CGM out of their halos, in contrast with other simulations. Finally, we predict that Civcolumn densities in the CGM of L_*galaxies are unlikely to depend on host galaxy SMBH mass. Our results show that the scatter in the low-mass end of the M_BH–σrelation may indicate how effective an SMBH is in the local redistribution of mass in its host galaxy. 

Abstract We examine the quenching characteristics of 328 isolated dwarf galaxies $\left({10}^8<M_{\mathrm{star}}/M_{\odot }<{10}^{10}\right)$within theRomulus25cosmological hydrodynamic simulation. Using mock-observation methods, we identify isolated dwarf galaxies with quenched star formation and make direct comparisons to the quenched fraction in the NASA Sloan Atlas (NSA). Similar to other cosmological simulations, we find a population of quenched, isolated dwarf galaxies belowM_star< 10⁹M_⊙not detected within the NSA. We find that the presence of massive black holes (MBHs) inRomulus25is largely responsible for the quenched, isolated dwarfs, while isolated dwarfs without an MBH are consistent with quiescent fractions observed in the field. Quenching occurs betweenz= 0.5–1, during which the available supply of star-forming gas is heated or evacuated by MBH feedback. Mergers or interactions seem to play little to no role in triggering the MBH feedback. At low stellar masses,M_star≲ 10^9.3M_⊙, quenching proceeds across several Gyr as the MBH slowly heats up gas in the central regions. At higher stellar masses,M_star≳ 10^9.3M_⊙, quenching occurs rapidly within 1 Gyr, with the MBH evacuating gas from the central few kpc of the galaxy and driving it to the outskirts of the halo. Our results indicate the possibility of substantial star formation suppression via MBH feedback within dwarf galaxies in the field. On the other hand, the apparent overquenching of dwarf galaxies due to MBH suggests that higher-resolution and/or better modeling is required for MBHs in dwarfs, and quenched fractions offer the opportunity to constrain current models. 

Abstract We study satellite counts and quenched fractions for satellites of Milky Way analogs inRomulus25, a large-volume cosmological hydrodynamic simulation. Depending on the definition of a Milky Way analog, we have between 66 and 97 Milky Way analogs inRomulus25, a 25 Mpc per-side uniform volume simulation. We use these analogs to quantify the effect of environment and host properties on satellite populations. We find that the number of satellites hosted by a Milky Way analog increases predominantly with host stellar mass, while environment, as measured by the distance to a Milky Way–mass or larger halo, may have a notable impact in high isolation. Similarly, we find that the satellite quenched fraction for our analogs also increases with host stellar mass, and potentially in higher-density environments. These results are robust for analogs within 3 Mpc of another Milky Way–mass or larger halo, the environmental parameter space where the bulk of our sample resides. We place these results in the context of observations through comparisons to the Exploration of Local VolumE Satellites and Satellites Around Galactic Analogs surveys. Our results are robust to changes in Milky Way analog selection criteria, including those that mimic observations. Finally, as our samples naturally include Milky Way–Andromeda pairs, we examine quenched fractions in pairs versus isolated systems. We find potential evidence, though not conclusive, that pairs, defined as being within 1 Mpc of another Milky Way–mass or larger halo, may have higher satellite quenched fractions. 

The interaction between supermassive black hole (SMBH) feedback and the circumgalactic medium (CGM) continues to be an open question in galaxy evolution. In our study, we use SPH simulations to explore the impact of SMBH feedback on galactic metal retention and the motion of metals and gas into and through the CGM of L ∗ galaxies. We examine 140 galaxies from the 25 Mpc cosmological volume, Romulus25, with stellar masses between 3 × 10 9 - 3 × 10 11 M ⊙ . We measure the fraction of metals remaining in the ISM and CGM of each galaxy, and calculate the expected mass of its SMBH based on the M−σ relation. The deviation of each SMBH from its expected mass, ΔMBH is compared to the potential of its host via σ . We find that SMBHs with accreted mass above the empirical M−σ relation are about 15\% more effective at removing metals from the ISM than under-massive SMBHs in star forming galaxies. Over-massive SMBHs suppress the overall star formation of their host galaxies and more effectively move metals from the ISM into the CGM. However, we see little evidence for the evacuation of gas from their halos, in contrast with other simulations. Finally, we predict that C IV column densities in the CGM of L ∗ galaxies may depend on host galaxy SMBH mass. Our results show that the scatter in the low mass end of M−σ relation may indicate how effective a SMBH is at the local redistribution of mass in its host galaxy. 

Abstract We explore the characteristics of actively accreting massive black holes (MBHs) within dwarf galaxies in the Romulus25cosmological hydrodynamic simulation. We examine the MBH occupation fraction, X-ray active fractions, and active galactic nucleus (AGN) scaling relations within dwarf galaxies of stellar mass 10⁸M_⊙<M_star< 10¹⁰M_⊙out to redshiftz= 2. In the local universe, the MBH occupation fraction is consistent with observed constraints, dropping below unity atM_star< 3 × 10¹⁰M_⊙,M₂₀₀< 3 × 10¹¹M_⊙. Local dwarf AGN in Romulus25follow observed scaling relations between AGN X-ray luminosity, stellar mass, and star formation rate, though they exhibit slightly higher active fractions and number densities than comparable X-ray observations. Sincez= 2, the MBH occupation fraction has decreased, the population of dwarf AGN has become overall less luminous, and as a result the overall number density of dwarf AGN has diminished. We predict the existence of a large population of MBHs in the local universe with low X-ray luminosities and high contamination from X-ray binaries and the hot interstellar medium that are undetectable by current X-ray surveys. These hidden MBHs make up 76% of all MBHs in local dwarf galaxies and include many MBHs that are undermassive relative to their host galaxy’s stellar mass. Their detection relies on not only greater instrument sensitivity but also better modeling of X-ray contaminants or multiwavelength surveys. Our results indicate that dwarf AGN were substantially more active in the past, despite having low luminosity today, and that future deep X-ray surveys may uncover many hidden MBHs in dwarf galaxies out to at leastz= 2. 

ABSTRACT We examine massive black hole (MBH) mergers and their associated gravitational wave signals from the large-volume cosmological simulation Astrid . Astrid includes galaxy formation and black hole models recently updated with an MBH seed population between 3 × 104h−1M⊙ and 3 × 105h−1M⊙ and a sub-grid dynamical friction (DF) model to follow the MBH dynamics down to 1.5 ckpc h−1. We calculate the initial eccentricities of MBH orbits directly from the simulation at kpc-scales, and find orbital eccentricities above 0.7 for most MBH pairs before the numerical merger. After approximating unresolved evolution on scales below $${\sim 200\, \text{pc}}$$, we find that the in-simulation DF on large scales accounts for more than half of the total orbital decay time ($$\sim 500\, \text{Myr}$$) due to DF. The binary hardening time is an order of magnitude longer than the DF time, especially for the seed-mass binaries (MBH < 2Mseed). As a result, only $$\lesssim 20{{\rm per \,cent}}$$ of seed MBH pairs merge at z > 3 after considering both unresolved DF evolution and binary hardening. These z > 3 seed-mass mergers are hosted in a biased population of galaxies with the highest stellar masses of $$\gt 10^9\, {\rm M}_\odot$$. With the higher initial eccentricity prediction from Astrid , we estimate an expected merger rate of 0.3−0.7 per year from the z > 3 MBH population. This is a factor of ∼7 higher than the prediction using the circular orbit assumption. The Laser Interferometer Space Antenna events are expected at a similar rate, and comprise $$\gtrsim 60\,{\rm{per\,cent}}$$ seed-seed mergers, $$\sim 30\,{\rm{per\,cent}}$$ involving only one seed-mass MBH, and $$\sim 10\,{\rm{per\,cent}}$$ mergers of non-seed MBHs. 

ABSTRACT In this work, we establish and test methods for implementing dynamical friction (DF) for massive black hole pairs that form in large volume cosmological hydrodynamical simulations that include galaxy formation and black hole growth. We verify our models and parameters both for individual black hole dynamics and for the black hole population in cosmological volumes. Using our model of DF from collisionless particles, black holes can effectively sink close to the galaxy centre, provided that the black hole’s dynamical mass is at least twice that of the lowest mass resolution particles in the simulation. Gas drag also plays a role in assisting the black holes’ orbital decay, but it is typically less effective than that from collisionless particles, especially after the first billion years of the black hole’s evolution. DF from gas becomes less than $$1{{\ \rm per\ cent}}$$ of DF from collisionless particles for BH masses >107 M⊙. Using our best DF model, we calculate the merger rate down to z = 1.1 using an Lbox = 35 Mpc h−1 simulation box. We predict ∼2 mergers per year for z > 1.1 peaking at z ∼ 2. These merger rates are within the range obtained in previous work using similar resolution hydrodynamical simulations. We show that the rate is enhanced by factor of ∼2 when DF is taken into account in the simulations compared to the no-DF run. This is due to $${\gt}40{{\ \rm per\ cent}}$$ more black holes reaching the centre of their host halo when DF is added. 

ABSTRACT We characterize the population of wandering black holes, defined as those physically offset from their halo centres, in the romulus cosmological simulations. Unlike most other currently available cosmological simulations, black holes are seeded based on local gas properties and are permitted to evolve dynamically without being fixed at halo centres. Tracking these black holes allows us to make robust predictions about the offset population. We find that the number of wandering black holes scales roughly linearly with the halo mass, such that we expect thousands of wandering black holes in galaxy cluster haloes. Locally, these wanderers account for around 10 per cent of the local black hole mass budget once seed masses are accounted for. Yet for higher redshifts ($$z$$ ≳ 4), wandering black holes both outweigh and outshine their central supermassive counterparts. Most wandering black holes, we find, remain close to the seed mass and originate from the centres of previously disrupted satellite galaxies. While most do not retain a resolved stellar counterpart, those that do are situated farther out at larger fractions of the virial radius. Wanderers with higher luminosities are preferentially at lower radius, more massive, and either closer to their host’s mid-planes or associated with a stellar overdensity. This analysis shows that our current census of supermassive black holes is incomplete and that a substantial population of off-centre wanderers likely exists.